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A chiral pion-nucleon amplitude supplemented by the HJS subthreshold coefficients is used to calculate the the long range part 
of the two-pion exchange nucleon-nucleon potential. In our expressions the HJS coefficients factor out, allowing a clear identi- 
fication of the origin of the various contributions. A discussion of the configuration space behaviour of the loop integrals that 
determine the potential is presented, with emphasis on cancellations associated with chiral symmetry. The profile function for 
the scalar-isoscalar component of the potential is produced and shown to disagree with those of several semi-phenomenological 
potentials. 
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I. INTRODUCTION 

Nowadays there are several semi-phenomenological NN 
potentials that may be considered as realistic because 
they provide good descriptions of cross sections, scatter- 
ing amplitudes and phase shifts. Nevertheless, it is pos- 
sible to notice important discrepancies when one com- 
pares directly their configuration space profile functions. 
Of course, this situation is consistent with the venera- 
ble inverse scattering problem, whereby there are always 
many potentials that can explain a given set of observ- 
ables. Therefore one must look elsewhere in order to 
assess the merits of the various possible models. 

In the case of NN interactions, there is a rather rich 
relationship between the potential and observables, in- 
volving several spin and isospin channels and different 
spatial regions. On the other hand, as all modern mod- 
els represent the long range interaction by means of the 
one pion exchange potential (OPEP), one must go to in- 
ner regions in order to unravel the discrepancies among 
the various approaches. 

Models vary widely in the way they treat the non- 
OPEP part of the interaction and, in the literature, one 
finds potentials constructed by means of dispersion rela- 
tions, field theory or just based on common sense guesses. 
In all cases, parameters are used which either reflect 
knowledge about other physical processes or are adjusted 
ad hoc. This leaves a wide space for personal whim and 
indicates the need of information with little model depen- 
dence about the inner part of the nuclear force. In the 
case of NN interactions, the complexity of the relevant 
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physical processes increases very rapidly as the internu- 
cleon distance decreases and hence the best process for 
yielding information with little model dependence is the 
tail of the two-pion exchange potential (TPEP). 

This problem has a long history. More than thirty 
years ago, Cottingham and Vinh Mau began a research 
program based on the idea that the TPEP is related 
to the pion-nucleon (ttN) amplitude ID. It lead to the 
construction of the Paris potential ^H, where the in- 
termediate part of the force is obtained from empirical 
ttN information treated by means of dispersion relations. 
This procedure minimizes the number of unnecessary hy- 
potheses and hence yields results which can be considered 
as model independent. Another important contribution 
was made by Brown and Durso Q who stressed, in the 
early seventies, that chiral symmetry has a main role in 
the description of the intermediate irN amplitude. 

In the last four years the interest in applications of 
chiral symmetry to nuclear problems was renewed and 
several authors have reconsidered the construction of 
the TPEP. At first, only systems containing pions 
and nucleons were studied, by means of non-linear la- 
grangians based on either PS or PV pion-nucleon cou- 
plings ^M0J|P]. Nowadays, the evaluation of this part 
of the potential in the framework of chiral symmetry has 
no important ambiguities and is quite well understood. 
This minimal TPEP fulfills the expectations from chi- 
ral symmetry and, in particular, reproduces at the nu- 
clear level the well known cancellations of the intermedi- 
ate ttN amplitude jHj|nj]. On the other hand, it fails to 
yield the qualitative features of the medium range scalar- 
isoscalar NN attraction p|p^]. This happens because a 
system containing just pions and nucleons cannot explain 
the experimental nN scattering data and one needs 
other degrees of freedom, especially those associated with 
the delta and the irN tr-term. The former possibility was 
considered by Ordonez, Ray and Van Kolck and 
shown to improve the predictive power of chiral cancel- 
lations but, in their work they did not regard closely the 
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experimental features of the intermediate ttN amplitude. 

Empirical information concerning the intermediate ttN 
process may be introduced into the TPEP in a model 
independent way, with the help of the Holder, Jacob 
and Strauss (HJS) subthreshold coefficients |T^,[l6|]. This 
kind of approach has already been extensively adopted 
in other problems. For instance, Tarrach and M.Ericson 
used it in their study of the relationship between nucleon 
polariz ability and nuclear Van der Waals forces JL7]|. In 
the case of three-body forces, it was employed in the con- 
struction of both model independent and model depen- 
dent two-pion exchange potentials |l^,[l^,|2(| . Using the 
same strategy, we have recently shown that the knowl- 
edge of the irN amplitude, constrained by both chiral 
symmetry and experimental information in the form of 
the HJS coefficients, provides a unambiguous and model 
independent determination of the long range part of the 
two-pion exchange NN potential [pl| . There we restricted 
ourselves to the general formulation of the problem and 
to the identification of the leading scalar-isoscalar poten- 
tial. In the present work we explore the numerical con- 
sequences of the expressions derived in that paper and 
compare them with some existing potentials. 

Our presentation is divided as follows: in Sec. II, we 
briefly summarize the derivation of the potential and rec- 
ollect the main formulae for the sake of self-consistency, 
leaving details to Appendices A, B, and C. In Sec. Ill we 
discuss the main features of the loop integrals that deter- 
mine the potential, emphasizing in the approximations 
associated with chiral symmetry. In Sec. IV we relate 
our theoretical expressions with those of other authors 
and in Sec. V, results are compared with existing phe- 
nomenological potentials. Finally in Sec. VI we present 
our conclusions. 



II. TWO-PION EXCHANGE POTENTIAL 

The construction of the ttttEP begins with the eval- 
uation of the amplitude for on-shell NN scattering due 
to the exchange of two pions. In order to avoid double 
counting, we must subtract the term corresponding to 
the iterated OPEP and, on the centre of mass of the NN 
system, the resulting amplitude is already the desired po- 
tential in momentum space. As it depends strongly on 
the momentum transferred A and little on the nucleon 
energy E, we denote it by T(A)u. 

In this work we are interested in the central, spin-spin, 
spin-orbit and tensor components of the configuration 
space potential, which may be written in terms of local 



*The final relativistic expression for the amplitude depends 
also on powers of z — p' + p, which yield "non-local" terms. 
We expand the amplitude in powers of z/m and keep just the 
first term, which gives the spin-orbit force. 



profile functions |2^]. They are related with the appro- 
priate amplitudes in momentum space by 



y ' \2mJ 4tt 



d J A 

(2^)3 
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In general, there are many processes that contribute 
to the TPEP. However, for large distances, the poten- 
tial is dominated by the low energy amplitude for 7rA 
scattering on each nucleon. When the external nucleons 
are on-shell, the amplitude for the process n a (k)N(p) — > 
TT b (k')N(p') is written as 



F = F+ 5 ab + F- i e bac t c 



(2) 



where 
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The functions and depend on the variables 

(p + p') ■ (k + k r ) 



t = (p — p') and v = 
or, alternatively, on 



4m 
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s = (p + k) 2 and u = [p — k') 2 



(5) 



When the pions are off-shell, they may also depend on 
k 2 and k' 2 . However, as discussed in Ref. off-shell 
pionic effects have short range and do not contribute to 
the asymptotic amplitudes. At low energies, A ± and B ± 
may be written as a sum of chiral contributions from the 
pure pion-nucleon sector, supplemented by a series in the 
variables v and t as follows: 
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In these expressions, the nucleon contributions were cal- 
culated using a non-linear pseudoscalar (PS) ttN cou- 
pling [ pT[ and, in writing A + , we have made explicit the 
factor (g 2 /m), associated with chiral symmetry. This 
amounts to just a redefinition of the usual a^~ , given in 
Refs. [^||,^6|. On the other hand, the use of a pseudo 
vector (PV) irN coupling would imply also a small redef- 
inition of b^Q. It is very important to note, however, that 
the values of the sub-amplitudes A ± and S ± are not at 
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FIG. 1. A) diagrams contributing to the low-energy ttN 
amplitude, where R represents the processes associated with 
the HJS coefficients; b) the two-pion exchange amplitude; c) 
contributions to the two-pion exchange amplitude from the 
purely pionic sector (top) and from processes involving the 
HJS coefficients (bottom). 



all influenced by this kind of choice and hence are com- 
pletely model independent. In the sequence, the terms 
in these expressions associated with the HJS coefficients 
will be denoted by and , the subscript R standing 
for "remainder" , as indicated in Fig. [l](A). 

The evaluation of the diagrams of Fig. |lj(B) yields the 
following general form for T 



T 
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2 J (2tt) 4 k 2 - ^ k' 2 - ^ 
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2t • 7" F-« F 



-(2) 
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where the are given in Eq. (||) and the factor \ ac- 
counts for the symmetry under the exchange of the in- 
termediate pions. The pion mass is represented by /i and 
the integration variable Q is defined as 

Q=\{k + k'). (11) 
In the sequence, we will also need the variables 

W=p 1 +p 2 =p' 1 +p' 2 , (12) 

A = k' - k = p[ - p x = p 2 - p' 2 , (13) 
Z = X - [(pi+pi)-(Pa +*£)]> (14) 
Vl = ±(W + Z ), (15) 



(16) 



The evaluation of the diagrams of Fig. |l|(C) produces 
1 f d 4 Q 
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where / and 
izations as 
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arc defined with non-rclativistic normal- 
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The integrand also depends implicitly on Q through the 
variables 



m 



ui — m = Q 
Vi = Q ■ Vi . 



Q-(W±z)-^ A 2 
■ Q.{W±z)-\ A 2 



(20) 

(21) 
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The integration is symmetric under the operation Q — > 
—Q and hence nucleon denominators involving s and u 
yield identical results. 

The evaluation of the potential in configuration space 
requires also an integration over t and the pole struc- 
ture of Eq. ( |Io[) implies that the leading contribution at 
very large distances comes from the region t « 4/i 2 p3[ , 
as it is well known. Therefore the form of our results 
in configuration space becomes more transparent when 
the contribution of the HJS coefficients is reorganized in 
terms of the dimensionless variable 



t 

v 



- 1 



(23) 



The amplitudes A^ and B^, associated with the HJS 
coefficients, are rewritten as 
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TABLE I. Values for the dimensionless coefficients of Eqs. ( fj^| - [§^ ) taken from Ref. J7^ / and re-stated by Eq. 



(m, n) 


(0,0) 


(0,1) 


(0,2) 


(1,0) 


(1,1) 


(2,0) 


mn 


3.676 ±0.138 


5.712 ± 0.096 


0.576 ± 0.048 


4.62 


-0.04 


1.2 ± 0.02 


0+ 


-2.98 ±0.10 


0.40 ± 0.04 


-0.16 


-0.68 ±0.06 


0.32 ± 0.04 


-0.31 ±0.02 




-10.566 ±0.212 


-1.976 ±0.144 


-0.240 ± 0.032 


1.222 ± 0.074 


0.208 ± 0.024 


-0.33 ±0.02 




9.730 ±0.172 


1.760 ± 0.104 


0.40 ± 0.032 


0.86 ±0.07 


0.22 ±0.02 


0.25 ±0.02 



.4 



(2m+l) 



(2m+l) 



(24) 
(25) 
(26) 
(27) 



In defining the coefficients (X^ n and we have introduced powers of // where appropriate so as to make them 

dimensionless. Their numerical values are given in Tab. |. 

Eq. (jlT]) can be naturally decomposed into a piece proportional to g 4 , which originates in the pure pion-nucleon 
sector and a remainder, labelled by R, as in Fig. 0(C). The former was discussed in detail in Refs. [p|Jl2[, where 
numerical expressions were produced, and will no longer be considered here. We concentrate on T^, which encompasses 
all the other dynamical effects. 

The potential in configuration space may be written as 

Vr = (V+ ± V+ + V+ + V+ A + V+ + V+ 6 + V+ + V+) + r« ■ (y R1 + V R2 + V m ± V m + V m + V m + V R7 + V RS ) 

(28) 



where the V R . are integrals of the form 
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and the gf are the polynomials in v/ fi and 6 given in Appendix A. Thus we obtain the following general result for 
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(35) 
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The expressions for V^ 3 , and are identical respectively to > Vr% an d V^f 7 when the very small differences 
between v\ and are neglected. In these results SB( m .n) an d ST( m .n) represent integrals of bubble (B) and triangle 
(T) diagrams, with m and n indicating the powers of (v/fj) and respectively, whose detailed form is presented in 
appendix B. There, we show that the integrals with one free Lorentz index are proportional to whereas those with 
two indices may be proportional to either V^V" or g^ u . Therefore we write for both bubble and triangle integrals 



cA 1 _ \rV qV 

D (m,n) — V i J (m,n)' 



Using the approximations described in Appendix B and the Dirac equation as in Eq.(Bl), we obtain 
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In configuration space, the spin-dependence of the potential is obtained by means of the non-relativistic results ]2S 



j(i)j(2) ^ x 
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where 
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An interesting feature of the partial contributions to 
the potential is that they are given by two sets of phe- 
nomenological parameters, the ttN coupling constant 
and the HJS coefficients, multiplying structure integrals. 
These integrals depend on just the pion and nucleon 
propagators and hence carry very little model depen- 
dence. Their main features are discussed in the next 
section. 




III. INTEGRALS AND CHIRAL SYMMETRY 



Our expressions for the TPEP, given by Eqs. (|42j-|49|), 
contain both bubble and triangle integrals, which depend 
on the indices m and n, associated respectively with the 
powers of {y/n) and 9, in the HJS expansion. The nu- 
merical evaluation of these integrals has shown that there 
is a marked hierarchy in their spatial behavior and that 
the functions with m = n = prevail at large distances. 

In order to provide a feeling for the distance scales of 
the various effects, in Figs. || and || we display the ratios 



qV 



/^T(0,0) 



for some values 



/Sb(o,o)] and 
of m and n, as functions of r. 

When considering these figures, it is useful to bear in 
mind that the (m, 0) and (0, n) curves convey different 
informations. The former series represents the average 
values of (y / '/i) m and is related to the behavior of the 
intermediate irN amplitude below threshold. For physi- 
cal ttN scattering, the variable v is always greater than 
fi, whereas in the present problem the average values of 
(yl\i) m are smaller than 1 for distances beyond 2.5 fm 
and tend to zero for very large values of r. This is the rea- 
son why the construction of the TPEP cannot be based 
on raw scattering data, but rather, requires the use of dis- 
persion relations in order to transform the ttN amplitude 
to the suitable kinematical region |2l| . One has, there- 
fore, a situation similar to the case of three-body forces, 
as discussed by Murphy and Coon |24| , which emphasizes 
the role of the HJS coefficients. 

Regarding the dependence of the integrals on the mo- 
mentum transferred, one notes that the intermediate ttN 
amplitude in the momentum space is already in the phys- 
ical t < region and does not require any extrapolations. 
On the other hand, when one goes to configuration space, 



FIG. 2. Asymptotic behavior of the bubble integrals 
5'_g(m,n)- The ratios S B ^ mi0 ^/S B ^ 0t0 ^ and SB(o,n)/^B(o,ff)! f or 
some values of m and nare indicated by solid and dashed lines 
respectively. One sees that the integral SWo.o) (unity line) is 
asymptotically dominant. 
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FIG. 3. Asymptotic behavior of the triangle integrals 
S T ( m ,n)- The ratios S T{mfi) /S T{0fi) and S T{0>n) / S T(ofl) , for 
some values of m and nare indicated by solid and dashed lines 
respectively. As in Fig. the integral for m — n — is 
asymptotically dominant. 
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FIG. 4. Structure of the leading contribution to the cen- 
tral potential, as given by Eq. pq. The continuous line 
represents the total effect, whereas the dashed, dotted, and 
dash-dotted lines correspond to the contributions proportional 

to (g 2 fi/m)a+ 00) 2S B , {g 2 n/m)a+ 00) 2St, and (a^ 00) ) S B re- 
spectively. 



the Fourier transform picks up values of the amplitude 
around the point t = Ay?. Thus, the r-space potential 
is not transparent as far as t is concerned and the co- 
herent physical picture only emerges when one uses it in 
the Schrodinger equation. This is a well known property 
which also applies to the OPEP. 

The fact that the integrals with m = n = dominate at 
large distances means that the main contribution to the 
isospin symmetric central potential comes from Eq. (|4^ ) 
and is given by 



Vri 



x 2 



V, 



R5 



fi 3 r 2 y + 



5' 



-8(0,0) 



+ ^T(0,0) 



^B(0, 0) > I 



(i)r(2) 



(55) 



The first term within curly brackets, proportional to 
g 2 , is produced by the triangle and bubble diagrams in 
Fig. [l](C)-bottom, containing nucleons on one side and 
HJS amplitudes on the other, whereas the second one is 
due to the last diagram of Fig. |l](C)-bottom. Inspecting 
Tab. g one learns that (g 2 ///m) /ckqo ~ 8, which suggests 
the first class of diagrams should dominate. On the other 
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-8(0,0) 



oV 

°T(0,0) 



hand, the first term is proportional to 

and, as discussed in appendix B, these two integrals have 
opposite signs and there is a partial cancellation between 
them. These features of the leading contribution are dis- 
played in Fig. ^, which shows that the first term is indeed 
dominant. 

The cancellation noticed in the leading contributions 
is not a coincidence. Instead, it represents a deep feature 
of the problem, which is due to chiral symmetry and also 
occurs in various other terms of the potential. 
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FIG. 5. Contributions of the box, crossed, and triangle di- 
agrams divided by that of the bubble, in the pure tyN sector, 
for the ratios of the pion over the nucleon mass equal to the 
experimental value fj,/m, to 0.1fj,/m, and to 0.01/x/m. 



In appendix C we have shown that the asymptotic form 
of 5^(0,0) is given by the analytic expression 



r*asymp 
,5 B(0,0) 



(4tt)« 
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(56) 



Its accuracy is 1% up to 1.2 fm. There, we also studied 
the form of the basic triangle integral St( ,o) an< ^ have 



demonstrated that S^qq? ~ 



0. 



S B(ofl) when (^/ m ) 
As the integrals with other values of m and n can be 
obtained from the leading ones, the same relationship 
holds for them as well. This explains why Figs. § and || 
are so similar. 

As we have discussed elsewhere [10 n|, important can- 
cellations due to chiral symmetry also occur in the pure 
ttN sector. In order to stress this point, we have evalu- 
ated the contributions of the diagrams in the top line of 
Fig. 0(C), denoted respectively by box (□), crossed (XI), 
triangle (A) (twice) and bubble (0), for three different 
values of the ratio y/m, namely 

Km) 10 \mJ 100 \mJ 

In Fig. H we display the ratios of the box, crossed and 
triangle contributions over the bubble result as functions 
of distance, where it is possible to notice two interesting 
features. The first is that these ratios tend to become flat 
as fi/m decreases. The other one is that as (n/m) — * 0, 
one obtains the following relations: □ = 0.50, 1X1= 0.5(), 
and A = — Q. Thus, for the amplitude in the pure ttN 
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sector, we have □+ IX +2A+0 = 0, a point also remarked 
by Friar and Coon H. This result, when combined with 
the previous discussion concerning the bottom part of 
Fig. |l](C), indicates that the two-pion exchange NN po- 
tential would vanish if chiral symmetry were exact, be- 
cause the same would happen with the intermediate ttN 
amplitude. So, all the physics associated with the tail of 
the intermediate range interaction is due to chiral sym- 
metry breaking. 

As a final comment, we would like to point out that 
in the evaluation of the TPEP there are two different 
hierarchies that can be used to simplify calculations. One 
of them concerns the HJS coefficients, which are more 
important for low powers on v and t. The other one 
is associated with the spatial behavior of the integrals 
as functions of m and n. The combined use of these 
hierarchies allow many terms to be discarded. 

IV. RELATED WORKS 

To our knowledge, only Ordonez, Ray and van Kolck 
have so far attempted to derive realistic nucleon-nucleon 
phenomenology in the framework of chiral symme- 
try |l4|,[lj|. The potential obtained by these authors 
is based on a very general effective Lagrangian, which 
is approximately invariant under chiral symmetry to a 
given order in non-relativistic momenta and pion mass. 
They considered explicitly the degrees of freedom associ- 
ated with pions, nucleons and deltas, whereas the effects 
of other interactions were incorporated into parameters 
arising from contact terms and higher order derivatives. 
In principle the free parameters in their effective La- 
grangian could be obtained from other physical-processes, 
but at present only some of them are knownlO. In their 
work these parameters were obtained by fitting deuteron 
properties and NN observables for j < 2 whereas loop in- 
tegrals were regularized by means of non-covariant Gaus- 
sian cutoffs of the order of the p meson mass. Thus they 
could show that the effective chiral Lagrangian approach 
is flexible enough for allowing the data to be reproduced 
with an appropriate choice of dynamical parameters and 
cutoffs. Comparing their approach to ours, one notes 
several important differences. For instance, we use di- 
mensional regularization which is well known to preserve 
the symmetries of the problem and our expressions are 
quite insensitive to short distance effects. In the work of 
Ordonez, Ray and van Kolck, on the other hand, "vari- 
ations in the cutoff are compensated to some exten t by 
a redefinition of the free parameters in the theory." p5| . 
Moreover, we use the HJS coefficients as input, which are 
determined by ttN scattering, and therefore our results 
yield predictions for interactions at large distances or, al- 



See Ref. j2a] for a comprehensive discussion of this point. 



ternatively, for j > 2. The test of these predictions will 
be presented elsewhere. 

Another point in the present work that deserves to be 
discussed concerns the subtraction of the iterated OPEP. 
In our calculation of the TPEP in the pure nucleonic sec- 
tor, we have supplemented the results derived by Lomon 
and Partovi [2^] for the pseudoscalar box and crossed 
box diagrams with bubble and triangle diagrams asso- 
ciated with chiral symmetry We have also shown 
that the use of a pseudovector coupling yields exactly 
the same results and hence that the potential does not 
depend on how the symmetry is implemented. However, 
the Partovi and Lomon amplitude include the subtrac- 
tion of the OPEP by means of the Blankenbecler-Sugar 
reduction of the relativistic equation and hence our re- 
sults are also affected by that procedure. This kind of 
choice should not influence measured quantities, since it 
amounts to just a selection of the conceptual basis to 
treat the problem [^6|. As discussed by Friar [^7| and 
more recently by Friar and Coon 0, the treatments of 
the iterated OPEP by Taketani, Machida and Ohnuma 
p8[ and by Brueckner and Watson [^9| differ by terms 
which are energy dependent. However, in our calcula- 
tion, energy dependent terms can be translated into the 
variable v and, in the previous section, we have shown 
that the TPEP at large distances is dominated by the 
region where v ss 0. Hence our results are not affected by 
the way the OPEP is defined. Another indication that 
confirms this fact comes from two recent studies dealing 
with the relative weights of the various TPEP contribu- 
tions to NN phase shifts, which have shown that the role 
of the iterated OPEP is very small for j > 2 

The last comment we would like to make in this sec- 
tion concerns the dynamical significance of the HJS coef- 
ficients. It has long been known that a tree model for the 
intermediate irN amplitude containing nucleons, deltas, 
rho mesons and an amplitude describing the cr-term can 
be made consistent with the experimental values of the 
HJS coefficients by means of a rather conservative choice 
of masses and coupling constants |l^,|24|,|3^|3f]|52| . In 
general, there are two advantages of employing such a 
model in a nuclear physics calculation. The first is that 
it allows one to go beyond the HJS coefficients, specially 
as far as the the pion off-shell behaviour of the ampli- 
tude is concerned. However, as we have discussed above, 
this kind of off-shell effects are related to short distance 
interactions and hence are not important for the asymp- 
totic TPEP. It is in this sense that we consider our 
results to be model independent. The second motivation 
for using a model is that it may provide a dynamical 
picture involving the various degrees of freedom of the 
problem and shed light into their relative importance. 
As we show in the next section, the leading contribution 
to the scalar-isoscalar potential comes from the coeffi- 
cient aQ = h(ciq + 4/z 2 aJ 1 + I6/j- 4 aJ 2 )- As expected, it is 
attractive and determined mostly by the irN sigma term 
and by the delta. The former yields a^ os = 1.8 whereas 
the latter is the outcome of a strong cancellation between 
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— __=-— 0.000 



2 3 



r (fm) 



FIG. 6. Structure of the central potential; the dot-dashed 
curve represents the leading contribution (Eq. fiplfy) whereas 
the dashed, big dotted and small dotted curves correspond to 
Eqs. (0), and Q) respectively; the solid line represent 

the full potential. 



pole and non-pole contributions ajj~ 0A = (26.5—25.2) p3[ . 
Thus the delta non-pole term plays a very important role 
in the interaction, and must be carefully considered in 
any model aiming at being realistic. 



V. RESULTS AND CONCLUSIONS 

In this work we have assumed that the TPEP is due 
to both pure pion-nucleon interactions and processes in- 
volving other degrees of freedom, as represented in the 
top and bottom lines of Fig. |l|(C). The former class of 
processes was evaluated and studied elsewhere [p 11 and 
hence we here concentrate on the latter. 

As discussed in Sec. Ill, the leading contribution to 
the potential at large distances is due to the intermedi- 
ate ttN amplitude around the point v = 0, t = 4/i 2 . In 
order to understand the role played by the other terms, in 
Fig. H we disclose the structure of the scalar- isoscalar po- 
tential, given by Eqs. ( p2[p^ ). There it is possible to see 
that Eq. (f42|), associated with the a+„ HJS coefficients, 
completely dominates the full potential. On the other 
hand, for moderate distances, there is a clear separation 
between the curves representing the leading contribution, 
given by Eq. ( |5q ) , and the total potential. This indicates 
that corrections associated with higher powers of v and 
t are important there, a feature that could have been 
anticipated from Figs. || and 

The total potential, obtained by adding the results of 
Refs. [||jl|] with those of this work, is given in Fig. |t], 
where it is possible to see that the contribution from the 
pure nucleon sector is rather small. This information, 
when combined with those contained in the preceding 
figures, allows one to conclude that the strength of the 
scalar-isoscalar attraction at large distances is due mostly 




r (fm) 

FIG. 7. Contributions for the total TPEP, represented by 
continuous line; the dashed line comes from, the pure nN sec- 
tor (Fig. ^(C)-top), whereas that associated with other degrees 
of freedom falls on top of the continuous line and cannot be 
distinguished from it. 



20 - ''\ 




r (fm) 



FIG. 8. Central components of various potentials: 
parametrized Paris |^/ (solid, P), Argonne vl4 J^ / (solid, 
A), dTRS (dashed, d), Bonn R|/ (dashed, B), and our 
full potential (solid, *). 



to diagrams involving the nucleon on one side and the 
remaining degrees of freedom on the other. 

In Fig. H we compare our results for the scalar-isoscalar 
interaction with the corresponding components of some 
potentials found in the literature: parametrized Paris [jj), 

Argonne vl4 §§, dTRS S and Bonn S The first 
thing that should be noted is that all curves but ours 
bend upwards close to the origin, indicating clearly that 
the validity of our results is restricted to large distances. 
Inspecting the medium and long distances regions, it is 
possible to see that every potential disagrees with all the 
others. On the other hand, this does not prevent the 
realistic potentials from reproducing experimental data, 
something that is possible because there is a compen- 
sation arising from the other discrepancies found in the 
short distance region. It is for this reason that the 
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2 3 4 5 6 7 

r (f m) 



FIG. 9. Ratio of the central components of some realistic potentials by our full result (solid,*): parametrized Paris jfy (solid, 
P), Argonne vl4 (solid, A), dTRS jj^ / (dashed, d), and Bonn |?|/ (dashed, B). 

accurate knowledge of the tail of the potential may yield indirect constraints over its short distance part. 

Finally, in Fig. Owe show the ratios of the realistic potentials by our full potential, where the discrepancies mentioned 
above appear again, in a different form. An interesting feature of this figure is that the realistic potentials come close 
together around 2 fm, suggesting that this region is important for reproduction of experimental data. Moreover, all of 
them show inflections there, indicating that the physics in this region goes beyond the exchange of two uncorrelated 
pions. In the long distance domain, the r dependence of the Argonne potential is not too different from ours, because 
it is based on a square OPEP form. 

In summary, in this work we have shown that the use of a chiral ttN amplitude, supplemented by experimental 
information, determines uniquely the long-distance features of the scalar-isoscalar component of the NN potential. 
As it is well known, the kinematical regions relevant to this problem are not directly accessible by experiment and 
hence empirical information has to be treated by means of dispersion relations before being used as input in the 
calculations of the force. From a purely mathematical point of view, our results are valid for r > 2.5 fm, since in this 
region one has v < // and the HJS coefficients may be safely employed. On the other hand, the determination of the 
dynamical validity of the results is much more difficult, since this requires a comparision with processes involving the 
mutual interaction of the exchanged pions, something that remains to be done in the framework of chiral symmetry. 

In general, a potential involves two complementary ingredients that deserve attention, namely geometry and dy- 
namics. In our calculation, the former is associated with standard bubble and triangle integrals, that determine 
unambiguously the profile functions in configuration space, whereas dynamics is incorporated into the problem by 
means of coupling constants and empirical coefficients. Geometry and dynamics decouple in our final expressions and 
hence they would remain valid even if changes in the values of the dynamical constants may occur in the future. In 
the case of Fig. |9j such a change would amount to just a modification of the vertical scale, with no appreciable effect 
on the discrepancies found with phenomenological potentials. 
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APPENDIX A: FUNCTIONS g? 



We present here the polynomials that enter Eq. (29). The groups of indices i = 1...4 and 5. ..8 refer, respectively to 
bubble and triangle diagrams. 
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The expressions for gf, g^r and g^ are obtained respectively from gf, g^ and gf by exchanging vx and v 2 



± „± 



APPENDIX B: INTEGRALS 



In this appendix we present the expressions for the integrals SB{ m ,n) and Sx( m ,n) that determine the potential 
given in Sec. II. In many cases, a considerable simplification of the results, with no loss of numerical accuracy, can 
be achieved due to the fact that one is interested in the asymptotic behaviour of the potential in configuration space. 
This allows one to ignore contact terms associated with delta-functions or, alternatively, constant terms in momentum 
space integrals. 

In bubble integrals the denominators involve just two pion propagators, whereas there is an extra nucleon propagator 
for triangles. In both cases, the integrands have the general form of a polynomial in the variables ^ = (® 1 VA , 
where Vi = -^{W ± z). In elastic pion-nucleon scattering at low energies, v = /i at threshold and hence [i is also a 
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natural unit for the Vi. In this problem, W « 2m, z p 
and hence Vi ~ 1. Moreover, using the mass-shell condi- 
tion for the external nucleons, we obtain 



(i) 



A 2 

4m 2 ' 

A 2 z 2 
4m 2 2m 2 



(Bl) 
(B2) 

(B3) 



These results show that the differences between v\ and v% 
are of the order of relativistic corrections and therefore 
may be neglected. 

In our expression for the potential, Lorentz tensors pro- 
portional to A's always appear contracted to either Vi or 
7 matrices. The use of the equations of motion imply 
in the vanishing of these products and hence we do not 
write them explicitly below. We also make use of the 
result A 2 = t = - A 2 . 

When going to configuration space, it is useful to use 
the following representation for the logarithm 
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1. Bubble integrals in momentum space: 

The basic bubble integral is 
d A Q 



1 B 



Q" 



(2tt)< [(Q_ f) 2 - M 2 ] [(Q + #) 2 -M 2 ]' 

(B5) 



The symmetry of the integrand makes all integrals with 
odd powers of Q to vanish. 

The simplest case corresponds to 



d A Q 



(2^) 4 [(Q_ | )2 _^] [ ( Q+A )2 _ m2 ] 

Using Feynman integration parameters we write 
d 4 Q 1 



Ib = 
where 



da 
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P =\ a -2 ) A ' 



M 2 =/ u 



2 -A 2 

4 



(B6) 

(B7) 

(B8) 
(B9) 



Using the technique of dimensional regularization de- 
scribed in 1 21 1, the integration over Q yields, after drop- 
ping the constant and divergent terms 
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Using Eq. (B4), we obtain 
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For the integral Ig , the same procedure yields 
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neglecting terms proportional to A M A y . 
Analogously, for 7g l/pA , we have 
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2. Triangle integrals in momentum space: 

The triangle integrals have the structure 

...... = f d A Q 2m^Sl...Sl 

T J (2tt) 4 [Q2 _ 2mV % ■ Q - ±A 2 ] 
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with 



P = ~[a-(1- a)(3] A - (1 - a)(l - /?)mV5, (B18) 

A 2 

M 2 = [a+ (1 — a) /3] p? + [1 - 2a - 2(1 - a)(3] — . (B19) 



Integrating over Q, we have 



where 
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Using the same procedure, and neglecting divergent terms, we get 
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3. Integrals in configuration space: 



The configuration space integrals are obtained by Fourier transforming the results given above multiplied by (— i) 
and by powers of the variable 6 = ^4^2 — l)- Recalling that t = — A 2 , we have the general structure 
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Neglecting contact terms, we obtain 
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where we use the short notation 
1 d 2 
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In general, the integrals that enter Eqs. (|3Q-p9D have 
at most two free Lorentz indices, since the other ones 
are contracted with powers of the vectors Vi . Therefore 
integrals with one free Lorentz index are proportional to 
and those with two indices are proportional to cither 
Vf V7 or g^ u , motivating the definitions of Eqs. 
In configuration space, the terms originating from the 
representation of the logarithm have the form 
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These integrals can be evaluated explicitly, and we 
have 
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For the bubble integrals this procedure yields the fol- 
lowing results 
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APPENDIX C: ANALYTICAL RESULTS FOR 
SOME INTEGRALS 

In this appendix we present analytic results for the 
asymptotic bubble and triangle integrals in configuration 
space needed in this work. 

The basic bubble integral is 
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For the triangle case, we have 
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where we have used a new variable s = (1 — /3)(1 — a)^. 
The function Mj. is given by eq.(B21) and can be rewrit- 
ten as 
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where y = vl + s 2 . Performing the y integration and 
comparing it with Eq. (CI), we find 
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Defining a new variable t such that 
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where x = fir. For large values of x, the integrand is very 
peaked around x rs and hence we expand the functions 
in front the exponential in a power series. Keeping the 
first three terms, we obtain our asymptotic expression 
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